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Abstract. - The Flory Huggins equation of state for monodisperse polymers can be turned 
into a density functional by adding a square gradient term, with a coefficient fixed by appeal 
to RPA (random phase approximation). We present instead a model nonlocal functional in 
which each polymer is replaced by a deterministic, penetrable particle of known shape. This 
reproduces the RPA and square gradient theories in the small deviation and/or weak gradient 
limits, and can readily be extended to polydisperse chains. The utility of the new functional is 
shown for the case of a polydisperse polymer solution at coexistence in a poor solvent. 



Introduction. - Recent work on polydisperse thermodynamics has shown the value of 
using, as thermodynamic coordinates, whatever linear combinations of densities actually occur 
in the excess free energy of the system. These linear combinations are 'generalized moments' 
(moments for short) of the complete set of densities, which are in principle infinite in number; 
however, in most models, particularly in polymeric systems, they form a much reduced set. 
For example, the Flory Huggins model of length-polydisperse homopolymers and solvent on a 
unit lattice has (setting kT = 1) F^ = (1 - <j>) ln(l - <f>) +X ( I>( 1 ~ where <f> = J dNNp(N) 
is the monomer density, or first moment of the distribution p(N) that specifies the density of 
chains of length TV. Thus there is only one moment in F ex . 

It was shown jl|-|3| how to construct a 'projected' free energy, expressed as a function of 
the required moments, whose predictions for cloud curves, shadow curves, and critical points 
coincide with those of the full model. It was also shown |^|,|| how, by systematically adding 
additional moments, the full phase behaviour (describing states of coexistence midway along 
a tie-ine as well as at its endpoints) could be recovered to any desired accuracy. 

An important extension of this approach is to address interface thermodynamic behaviour, 
for example to predict the interfacial tensions between coexisting phases. (In doing this, we can 
suppose that the composition of such phases was calculated already, e.g., by the methods just 
outlined.) This requires a model, not merely of the equation of state, but of the free energy 
functional. In some cases (such as polydisperse hard spheres) accepted functionals include 
ones where the excess free energy depend only on a few nonlocal 'weighted densities'. The 
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minimization equations then close in the subspace of these weighted densities Q , simplifying 
practical calculation. So far, there is no such treatment of polydisperse polymers, despite the 
simplicity of the excess free energy. This is partly because polymers are extended objects: the 
main difficulty lies in generalizing the ideal (entropic) term to inhomogeneous states. 

RPA and Square Gradient Models. - For concentrated monodisperse (homo)polymers, 
there are two simple methods for generating density functionals. The first is to expand in the 
random phase approximation (RPA) for small deviations about a uniform reference state. In 
Fourier space this gives to quadratic order 

AF=i5>(g)| 2 (So(q)- 1 +FM) (1) 

So(q) = <W<z) = p a N 2 g D (q 2 R 2 J ) (2) 

Here vector notation for the wavevector is suppressed and F" x is evaluated in a uniform 
state, of monomer density 0o = poN. The unperturbed scattering function So(q) describes 
an ideal gas of Gaussian chains (each with structure factor so(q)) at number density p. In 
cq.^j, go(x) — 2(exp(~x) — 1 + x)/x 2 is the Dcbye function and R g the chain gyration radius 
(on a unit lattice, (iV/6) 1 / 2 ) fj. Within RPA, one can continue to higher order, but with 
increasingly complicated results fj] . A disadvantage of eqjl] is that it does not allow study of 
(say) the interface between phases at coexistence, whose composition difference is finite. 

The second (square gradient) method allows for such finite deviations but assumes 
instead weak gradients, asserting that (again suppressing vector notation) 

AF = J [F(<j)(r))+K((t>)\\/<j)\ 2 ] dr (3) 

where F(<f>) = F id + F ex ((f>) and F id = (cf> / N)[ln(4>/ N) - 1] is the ideal contribution to the 
Helmholtz free energy of a homogenous phase. The coefficient K((f>) is a priori unknown, but 
can be identified by comparing with eq.[j]at order q 2 . This gives (on a unit lattice) K~ x = 360. 
This choice for K ensures that RPA is recovered when <f) variations are small as well as slowly 
varying. For strong segregation (when neither applies) some authors (but not others Q) 
propose matching instead the coefficient of |V0| 2 against the high-g limit of So(q)^ 1 which is 
also quadratic |9|; this gives K^ 1 = 240. 

An Alternative Functional. - The implementation of either of these methods for poly- 
disperse chains is less simple than it looks (see and below). We offer instead a density 
functional for dense polymer systems which reduces to the RPA and square gradient func- 
tionals in the appropriate limits, but whose extension to polydispersity is both transparent 
and simple to implement. In effect, we replace each chain by a nondeformable 'effective parti- 
cle' with extremely non-pairwise-additive interactions. (This last feature contrasts with, and 
complements, a related recent approach to dilute polymers p"l|].) Our functional interpolates, 
in an illuminating way, the physics of RPA and of the square gradient limits. 

For a single polymeric species in a monomeric solvent, we propose the model functional 

AF = J {p{r)[\np(r)-l]+F ex {4>{r))} dr (4) 

p(q) = (Nsoiq))- 1 / 2 <t>(q) (5) 

Expanding this functional in small deviations about a uniform reference state recovers imme- 
diately eq.|l], whereas expanding it in gradients yields eq.[| The latter expansion is achieved by 
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writing p(r) = (1/AQfl — ( TV/36) V 2 )0(r) (this is, in real space, the low q expansion of eq||), 
substituting into eq.^ and integrating by parts. The result K~ x = 360 is recovered without 
further assumption. These facts are sufficient to ensure that our functional encompasses, for 
monodisperse chains, the standard (leading order) RPA and square-gradient functionals. 

Interpretation. - Eq.^ has a simple physical interpretation, as follows. It describes a fluid 
of identical, interpenetrable, spherical particles of fixed structure. The density of the centres of 
these particles is p(r) and the ideal part of the free energy functional is simply p(r) [In p(r) — 1] 
(as applies for any simple fluid). Viewing each particle as a nondeformable cloud of monomers, 
it follows that the monomer density (f> is constructed from p by a deterministic, invertible real- 
space convolution, or Fourier product (eq.^|): 4>{q) — w(q)p(q). The kernel w is, remarkably, 
not the Debye function for a Gaussian chain, but proportional to its square root: w(q) = 

1/2 

Ngp . This choice ensures that the monomer-monomer correlator (in Fourier space, So(q)) 
in an ideal gas of our penetrable spheres is identical to that of an ideal gas of Gaussian chains. 
The kernel w itself is not the monomer-monomer correlator of an effective particle, but (in 
real space) the conditional probability of finding a monomer at r given that the centre of 
our effective particle is at the origin (f2) . Each effective particle is a spherically pre-averaged 
polymer in the sense that it has the same monomer-monomer correlator; but the density 
profile about its centre is not that about any point on a real chain (nor the centre of mass) . 
Indeed, since a polymer is stochastic, one cannot have both things at once; we give priority 
to So(q)- Finally (as in both the RPA and the square gradient models) F ex is taken as a 
local functional of the monomer density <f>. This reflects the fact that this term arises from 
enthalpic interactions unrelated to chain connectivity, whose scale is local (monomeric) Q. 

Eq.[I| illustrates the central issue in constructing density functionals for extended objects. 
Within Flory Huggins, the entropy in any homogeneous state is that of an ideal gas of chains. 
But a chain is not pointlike, and sees variations in density - so, which density should we use? 
The 'effective particle' description gives a definite (though perhaps not unique) answer. If 
eqji] is written down only in terms of </>, the entropy term is nonlocal Q . If written only in 
terms of p, then the enthalpic term is nonlocal instead. This 'uncertainty principle' will hold 
for any theory of extended objects with local interactions. Both p and <f> are needed, and our 
approximation lies in adopting a fixed relation between them, instead of the actual one which 
for polymers is (a) stochastic and (b) subject to perturbation by strong density variations. 

Our effective particles cannot adapt to their surroundings by deformation, only by dis- 
placement; this limits the formal range of validity for eq.|| mainly to cases where either RPA 
or the square gradient theory could have been validly used instead. On the other hand, by 
adopting a square gradient form also in strong segregation ||, the nondeformability of our 
particles could be relaxed, which might better approximate the behaviour of real polymers (to 
the same extent that strong segregation would do). In what follows we consider only weak 



,V L 

and intermediate (i? g V ~ 1) segregation where eq.[|is at least competitive. 

Polydispersity. - The present approach comes into its own in the presence of polydis- 
persity. The underlying 'effective particle' interpretation makes the required extension both 
obvious and tractable. To model lcngth-polydispersc polymers at number density p(N) (where 
N is now treated as a continuous variable) and with an excess free energy that depends (as 
before) only on the local monomer density <f>{r), we must choose 



AF= I dr{ I p(N,r)[lnp(N,r) - l]dN + F ex (4>(r))j (6) 

cf>(N, r) dN (7) 
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<j)(N,r) = J w(N,r-r')p(N,r')dr' (8) 

where, in Fourier space, w(N,q) = Ngu(q 2 N/6) 1 ' 2 , as before. The convolution between 
4>(N, r) (the monomer density at r from objects with N monomers) and p(N, r) (the density 
of centres of these objects) remains invertible; thus eq.|^ directly generalizes eq.|B|. 

Eq.^ allows practical numerics, because it is in the class where the excess free energy 
depends on only a few 'weighted densities' or nonlocal linear combinations of the particle 
densities p(N,r). In such cases, it is possible (once the bulk chemical potentials are known 
- equivalent to specifying the bulk composition po(N) of any one of the phases present) to 
exactly obtain autonomous (closed) equations among these weighted densities g. The size 
distribution at any point in the density profile can be calculated afterwards. 

In the present case there is only one weighted density (0) with the autonomous equation: 



w(N, r 



Po(iV)exp (- J w{N,r" -r')AF' ex [<j>{r"))dr"\ 



dNdr (9) 



where AF e ' 2 .(0) = F' ex {4>) — Fg X ((f>o) in which the prime denotes a derivative; 4>o is the bulk 
monomer density in the chosen reference phase. Note that AF^ x (<j)) is a function only of cj>, 
hence the closure at this level; however the integral over N means that eq.|| is not (to our 
knowledge) the Euler-Lagrange equation from any functional of <p(r) alone. 

Example. An example follows of a practical application of the method, generating 
interfacial profiles and properties for coexisting phases of polydisperse homopolymer in a 
solvent . We have considered, for simplicity, the case where the parent phase has a uniform 
length distribution (with mean length (N)q = 1000 and standard deviation 346 in the parent). 
We study coexistence between a cloud phase and its shadow, solving eq.^; the cloud phase 
is then the parent or reference po(N) and we find the interface between it and the (unique) 
shadow phase at a given x parameter. (By using an 'effective parent' the same method extends 
to any point along a tie line; see ||.) In fig.|l|a, we show the relevant coexistence curves, which 
display a dense polymer solution (cloud) coexisting with a dilute one (shadow). The length 
distribution in the shadow differs from the (flat) distribution in the cloud, as shown in figjljb; 
for our flat parent the size distribution in the shadow is exponential. The shorter chains are 
preferentially segregated into the shadow phase, which, in addition, has fewer chains overall. 

We have solved numerically eq|9j which gives us access to both structural and thermody- 
namic quantities for the interface between coexisting phases. For numerical convenience, the 
Debye function was replaced by g(x) = [1 + x 2 /2}~ 1 as is standard practice fllj. The first 
moment density (or total monomer concentration) exhibits smooth, monotonic profiles. From 
these, we can define an interfacial width, which shows a smooth decrease when increasing the 
interaction strength, \, However, at the same time we find a nonmonotonic distribution of 
monomer concentrations for individual species when we move from the vicinity of the critical 
point. In fig. ||a, we show the monomer concentration profiles contributed by various different 
chain lengths, centered at the position where <f>(r) takes its mean value. The shorter polymers 
show a maximum in the monomer concentration at the interface: they act rather like a sur- 
factant. For the same selection of species, we show in fig.||b profiles for the relative amounts 
of each species, c(N, x) = 4>{N, x)/4>{x). The strong depletion of the long chains in the dilute 
phase induces a nonmonotonic, even oscillatory, behaviour in this quantity. Closer to the 
critical point, the profiles for the monomer and relative concentrations change monotonously 
from the dense to the dilute phase. 

In figj^a we show the monomeric adsorption as a function of chain length for a fixed parent 
at various <fi values. The adsorbed amounts are defined with respect to the Gibbs dividing 
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surface of the overall monomer concentration. The smaller chains are adsorbed preferentially 
(consistent with fig.^a), and the species that has the maximum adsorption corresponds to a 
smaller chain for stronger demixing. The latter can be attributed to an increased skewness 
of the shape distribution in the dilute phase Jlij . Finally, in fig.|]b we compare the surface 
tension of the polydisperse fluid with an equivalent monodisperse one. It shows the surfactant 
effect of the short chains - the surface tension of the polydisperse mixture increases more 
slowly than its monodisperse counterpart - though the effect is surprisingly weak. A more 
detail analysis of these and the other results will be pursued elsewhere [Q . 

RPA and Square Gradient Limits. - Having shown its utility by this example, we re- 
turn to the general properties of our functional. Expanding in small deviations (RPA) now 
allows two separate descriptions, according to whether the densities of individual species are 
eliminated or retained. The simplest RPA is to keep eqjl] but to replace eq.|| with the ideal 
structure factor of the polydisperse noninteracting system (see, e.g. Jl^]). This is surely cor- 
rect, but allows no calculation of partial structure factors for differently labelled species and 
cannot be extended to dynamics [jUj , or even interfaces. 

The full RPA limit of our functional is instead 

AF =\Y, f IWril'fSo^r'+^W) dN (10) 

q=£0 J 

This coincides with the second of two distinct polydisperse RPAs suggested by Warren |l(J ; for 
our effective particles (at least) this is rigorous, and his first suggestion, which treats p{N, r) 
and <fi(N,r) as though independent, incorrect. 

Turning to the square gradient theory, we obtain from our functional (after expansion and 
integration by parts) 



AF = dr 



[—) r {—) - 



1 (V0(iV,r))H dN + F ex (cj)(r)) 



36</)(N, r) 



(11) 

as used, for example, in |L6|]. This is certainly the proper extension of eq.||, although the fact 
that the resulting equations can be closed in 4>{r) has, we think, been overlooked previously. 
In fact, at square gradient level eq|9| is converted from an integral equation into a second 
order differential equation for <p(r). However, this is so nonlinear that it offers little numerical 
advantage over eq.|^ and we found it just as easy to use the full functional. Analytically, the 
expansion in gradients does allow some limiting trends to be extracted |l4| ; a square-gradient 
truncation might also be used to mimic chain deformability in the strong segregation limit || . 

Extensions. - No assumption about F ex was made other that it was a local function of 
the overall monomer density <f>(r). The work extends straightforwardly to cases where F ex 
depends on more than one weighted density of p(N,r), for example the case of 'chemical 
polydispersity' (a model of random copolymers) in which chains have the same length but 
varying interaction parameter \ IMI - (This requires two weighted densities [gj.) Polydisperse 
polymer mixtures of several chemical types can likewise be included; the solvent can be made 
polymeric or removed if desired (by imposing (f> = 1 as a constraint). 

The work of Ref. || shows that excess free energies of up to four moments, at least, can 
be treated comfortably, which is probably more than required for most interfacial problems 
involving dense polymers. 

Conclusion. - We have introduced, within a Flory-Huggins framework, a model density 
functional for polydisperse polymers. This should be valid quantitatively in weak segregation 
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and qualitatively for moderate segregation. By representing chains as a system of spherically 
symmetric, nondeformable but interpenetrating particles, the functional is able to reproduce 
the known limiting forms for RPA and square gradient theories in the monodisperse limit. 
To study the interface thermodynamics of coexisting phases, the closure that we obtain at 
first-moment level (eq.^ ensures numerical tractability for our model functional, eq.^. This 
closure is inherited by the square gradient expansion of our functional, which has the standard 
polydisperse form (eq.|lT|); but in this context the existence of a closed differential equation 
for the monomer density has perhaps been overlooked. The practical utility of eq.^ was 
demonstrated by applying it to find interfacial tensions and density profiles for a polydisperse 
dense polymer solution at coexistence in a poor solvent. 

* * * 
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Fig. 1 - (a) Cloud (circles) and shadow (filled squares) curves for a uniform chain length distribution. 
The continuous curve is the spinodal. (b) Composition in the shadow phase for various points on the 
cloud curve (identified by <f> in the cloud phase, whose dependence on x is shown in (a)). 
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Fig. 2 - Monomeric density profiles for various of species for a flat parent at its cloud point, with 
{4>}o = 0.16 and x/Xc = 1-07. (a) Monomer concentration profiles, each normalized by its value in 
the cloud phase, (b) Relative concentration profiles (for the same species). 
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Fig. 3 - (a) Monomeric adsorption F, as a function of chain length, for the same parent as in the 
previous figures, at various points on the cloud curve, (b) Surface tension for the flat parent and an 
equivalent monodisperse solution, i.e. with the same < N >o. 



